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(The BaBAR Cohaboration) 

We measure the mass difference, Amo, between the Z)*(2010)''' and the Z)° and the natural line 
width, r, of the transition _D*(2010)'*" -^ D"tt~^ . The data were recorded with the BABAR detector 
at center-of-mass energies at and near the T{4S) resonance, and correspond to an integrated lumi- 
nosity of approximately 477 fb~^. The _D'^ is reconstructed in the decay modes D" — >■ K~7v^ and 
D° -^ A'~7r+7r"7r+. For the decay mode D° -s- A'"7r+ we obtain F = (83.4 ± 1.7 ± 1.5) keV 
and Amo = (145 425.6 ± 0.6 ± 1.8) keV, where the quoted errors are statistical and system- 
atic, respectively. For the D^ -^ K^tv^tt^tv^ mode we obtain F — (83.2 ± 1.5 ± 2.6) keV and 
Amo = (145 426.6 ± 0.5 ± 2.0) keV. The combined measurements yield F = (83.3 ± 1.3 ± 1.4) keV 
and Attio = (145 425.8 ± 0.5 ± 1.8) keV; the width is a factor of approximately 12 times more precise 
than the best previous value, while the mass difference is a factor of approximately 6 times more 
precise. 

PACS numbers: 13.20.Fc, 13.25Ft, 14.40.Lb, 12.38.Gc, 12.38.Qk, 12.39.Ki, 12.39. Pn 



INTRODUCTION 

The D*(2010)+ (£>*+) line width provides a window 
into a nonpcrturbativc regime of strong physics where 
the charm quark is the heaviest meson constituent |l|- 
|3|. The line width provides an experimental check of 
models of the D meson spectrum, and is related to the 
strong coupling of the Z?*"*" to the Dn system, gwDi^- 
In the heavy-quark limit, which is not necessarily a good 
approximation for the charm quark [J| , this coupling can 
be related to the universal coupling of heavy mesons to 
a pion, g. There is no direct experimental window on 
the corresponding coupling in the B system, gs'Sin since 
there is no phase space for the decay B* -> Btt. However, 
the D and B systems can be related through g, which 
allows the calculation of gs-'Bw- The B*Bt: coupling is 
needed for a model- independent extraction of |Vu&| |5|,|6| 
and is presently one of the largest contributions to the 
theoretical uncertainty on \Vub\ [Z|. 

We study the _D*+ -^ L)"7r+ transition using the D" -^ 



K^TT^ and D" -^ A'^7r+7r^7r+ decay modes to measure 
the values of the D*+ line width, F, and the difference 
between the D*^ and D" masses, Atoq [8|. The best 
prior measurement of the width is F = (96 ± 4 ± 22) keV 
by the CLEO collaboration where the uncertainties are 
statistical and systematic, respectively [9|. That mea- 
surement is based on a data sample corresponding to 
an integrated luminosity of 9fb^ and reconstructed 
D'^ — > K^n^ decays. In the present analysis, we have a 
data sample that is approximately 50 times larger. This 
allows us to apply tight selection criteria to reduce back- 
ground, and to investigate sources of systematic uncer- 
tainty with high precision. 

The signal is described by a relativistic Breit-Wigner 
(RBW) function defined by 



dF(m) 
dm 



(m2 
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where Fd'Dtt is the partial width to D^tt^, m is the 
D^TT^ invariant mass, ttiq is the invariant mass at the 



pole, and rTotai('TT-) is the total £)*+ decay width. The 
partial width is defined by 
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where J^^^ (p) — y/T+r^^ is the Blatt-Weisskopf form 
factor for a vector particle with radius parameter r and 
daughter momentum p, and the subs crip t zero denotes 
a quantity measured at the pole [10|, lll| . The value of 
the radius is unknown, but for the charm sector it is 
expected to be '-^ 1 GeV~^ [l^- We use the value r ~ 
1.6GeV~^ from Ref. [13[ and vary this value as part of 
our investigation of systematic uncertainties. 

Key to this analysis is that the tails of the RBW are 
much longer than the almost Gaussian resolution func- 
tion, which allows us to use information far from the 
central signal region to measure the width. 

This paper is organized as follows. Section n discusses 
the BaBAR detector and the data used in this analysis, 
and Section „ describes the event selection. Section n dis- 
cusses a correction to the detector material model and 
magnetic field map. Section „ details the fit strategy. Sec- 
tion n discusses and quantifies the sources of systematic 
uncertainty, and Section „ describes how the results for 
the two D° decay modes are combined to obtain the final 
results. Finally, the results are summarized in Section „ 



THE BABAR DETECTOR AND DATA 

This analysis is based on a data sample correspond- 
ing to an integrated luminosity of approximately 477 fb^ 
recorded at and 40 MeV below the T (45) resonance by 
the BaBAR detector at the PEP-II asymmetric energy 
e+e" collider. The BABAR detector is described in detail 



elsewhere 1^ , so we summarize only the relevant compo- 
nents below. Charged particles are measured with a com- 
bination of a 40-layer cylindrical drift chamber (DCH) 
and a 5-layer double-sided silicon vertex tracker (SVT), 
both operating within the 1.5-T magnetic field of a su- 
perconducting solenoid. Information from a ring-imaging 
Cherenkov detector is combined with specific ionization 
(dE/dx) measurements from the SVT and DCH to iden- 
tify charged kaon and pion candidates. Electrons are 
identified, and photons measured, with a CsI(Tl) elec- 
tromagnetic calorimeter. The return yoke of the super- 
conducting coil is instrumented with tracking chambers 
for the identification of muons. 



EVENT SELECTION 

We reconstruct continuum-produced Z?*+ -^ D^T^t ^'^~ 
cays in the two Cabibbo- favored channels D^ -^ K^tt^ 



and Z?" -^ K~Tr~^TT~TT^ . The pion from the D*~^ decay is 
called the "slow pion" (denoted 7r+) because of the lim- 
ited phase space available. The selection criteria for the 
individual D" channels are detailed below; however, both 
modes have the same D*~^ requirements. The selection 
criteria were chosen to enhance the signal-to-background 
ratio (S/B) to increase the sensitivity to the long RBW 
tails in the Am distribution; we have not optimized the 
criteria for statistical significance. Because this analysis 
depends on the RBW tails, we pay particular attention 
to how the selection criteria affect the tail regions. 

The entire decay chain is fit using a kinematic fitter 
with geometric constraints at each vertex, and the addi- 
tional constraint that the D*^ candidates originate from 
the interaction region. The confidence level of the x^ 
of the fit to the entire decay chain is required to be at 
least 0.1%. If there is more than one D*+ candidate in 
the event, we choose the one with the highest confidence 
level. The mass difference between the D*^ and Z?" is 
required to satisfy Am < 0.17 GeV. A large amount of 
the combinatorial and B meson decay background is re- 
moved by requiring p*(D*^) > 3.6 GeV, where p* is the 
momentum measured in the e~^e~ center-of-mass frame 
for the event. Selecting Cabibbo-favored decay modes 
avoids a large amount of combinatorial background mim- 
icking signal from doubly Cabibbo-suppressed decays. 

The experimental uncertainty in the measured invari- 
ant mass difference is dominated by the uncertainty on 
the measured momentum of the slow pion. Therefore, 
we choose criteria to select well-measured pious. For im- 
proved tracking resolution we require that the 7r+ tracks 
have at least 12 measurements in the DCH and have at 
least 6 SVT measurements with at least 2 in the first 
three layers. We reduce backgrounds from other species 
of tracks in our slow pion sample by requiring that the 
dE/dx values reported by the SVT and DCH be consis- 
tent with the pion hypothesis. 

The Dalitz decay n^ — >■ ^e'^e~ results in a problematic 
background where we misidentify an electron as a tt^ . We 
eliminate such candidates by reconstructing a candidate 
e+e" pair and combining it with a 7. If the e+e" vertex 
is within the SVT volume and the invariant mass is in the 
range 115 McV < m{'ye~^e~) < 155 MeV, then the event 
is rejected. Photon conversions in the detector material 
are another source of background where electrons can be 
misidcntificd as slow pious. To identify such conversions 
we first create a candidate e+e" pair using the slow pion 
candidate and an identified electron track from the same 
event and perform a least-squares fit with a geometric 
constraint. The event is rejected if the invariant mass of 
the putative pair is less than 60 MeV and the constrained 
vertex position is within the SVT tracking volume. 

We identified additional criteria to remove candidates 
in kinematic regions where the Monte Carlo (MC) sim- 
ulation poorly models the data. The MC is a cock- 
tail of qq and £~^£~ sources where q = u,d,s,c,b and 



i = e,iJ,, T. The simulation does not accurately replicate 
the momentum distributions observed in data at very 
high and low !)*+ momentum values, so we require that 
3.6 GeV < p*{D*+) < 4.3 GeV and that the laboratory 
momentum of the slow pion be at least 150 MeV. In an 
independent sample of Kg — )■ tt^tt^ decays, the recon- 
structed Kg mass is observed to vary as a function of 
the polar angle of the Kg momentum measured in the 
laboratory frame with respect to the electron beam axis. 
We define the acceptance angle to reject events where 
any of the daughter tracks of the D*+ has cos 6* > 0.89 
to exclude the very- forward region of the detector. This 
criterion reduces the final data samples by approximately 
10%. 



For both £>" 



K 7r+ and D'^ -^ K tt+tt tt"*", 



apply particle identification (PID) requirements to the 
K and tt candidate tracks. To select candidates with 
better tracking resolution, and consequently improve the 
resolution of the reconstructed masses, we require that 
D'^ daughter tracks have at least 21 measurements in the 
DCH and satisfy the same SVT measurement require- 
ments as the slow pion track described before. The re- 
constructed D" mass must be within the range 1.86 GeV 
to 1.87 GeV, and the confidence level of the x^ oi the fit 
of the D'^ with kinematic and geometric constraints must 
be at least 0.5%. 

The background level in the I?° -^ K^tt^tt^tt^ mode 
is much higher than that in D*^ — >■ K~Tr^, and so we 
require D° daughter charged tracks to satisfy stricter 
PID requirements. The higher background arises be- 
cause the D° mass is on the tail of two-body K~Tr^ 
invariant mass distribution expected in a longitudinal 
phase space model, however it is near the peak of the 4- 
body if ~7r"'"7r~7r"'" invariant mass distribution |15| . Com- 
pounding this phase space effect is the greater probabil- 
ity of wrongly assigning a track while still mimicking a 
correctly reconstructed 1?°. 

The initial fit to the D^ — > K^tt^tt^tt^ validation sig- 
nal MC sample had a bias in the measured value of the 
_D*+ width. An extensive comparison revealed that the 
bias originated from regions of phase space that the MC 
generator populated more frequently than the data. Ev- 
idently, there are amplitudes that suppress these struc- 
tures in the data, that arc neither known nor included 
in the MC generator. We avoid the regions where the 
MC disagrees with the data by rejecting a candidate if 
either (m^ (7r+7r+) < -1.17 m^ (7r-7r+) -H 0.46 GeV^) or 
(to2 (7r-7r+) < 0.35 GeV^ and m^ {K-t:+) < 0.6 GeV^). 
This veto is applied for each 7r+ daughter of the D" can- 
didate. 

There is an additional source of background that must 
be taken into account for the K~'k~^t:~t:^ channel that 
is negligible for the K~Tr~^ channel. In a small fraction 
of events (< 1%) we mistakenly exchange the slow pion 
from £)*+ decay with one of the same-sign _D° daughter 
pions. This shifts the reconstructed mass values and in- 



troduces a 0(0.1 keV) bias on the width obtained from 
the fits to the validation signal MC sample. To veto these 
events we recalculate the invariant mass values after in- 
tentionally switching the same-sign pions, and create the 
variable Am' = m(K~Tr~^7T~'!T^'!Tf) — m{K~'K^'K~'iTf). 
There are two pions from the D'^ decay with the same 
charge as the slow pion, so there are two values of Am! 
to consider. In this procedure the correctly reconstructed 
events are moved away from the signal region, while 
events with this mis-reconstruction are shifted into the 
signal region. We reject events with Am! < 0.1665 GeV. 
Using fits to the validation signal MC sample we found 
that by using this criterion rather than using a narrow 
range, we are able to remove approximately 80% of the 
misreconstructed events and this removes the bias of the 
fitted value of the width. 



MATERIAL MODELING 

In the initial fits to data, we observed a very strong 
dependence of the RBW pole position on the slow pion 
momentum. This dependence is not replicated in the 
MC, and originates in the magnetic field map and in the 
modeling of the material of the beam pipe and the SVT. 
Previous BaBAR analyses have observed such effects. For 
example, in the measurement of the 4+ mass [l6| the 
material model of the SVT was altered in an attempt 
to correct for the energy loss and the under-represented 
small-angle multiple scattering (due to nuclear Coulomb 
scattering). However, the momentum dependence could 
be removed only by adding an unphysical amount of ma- 
terial to the SVT. The corrections described here avoid 
altering the underlying detector model and focus on ad- 
justing track momenta after reconstruction. 

We use a sample of Kg — s- tt+tt" candidates from 
£)*+ — 7> Z?°7r+ decay, where we reconstruct D° — )■ 
KgiT'^TT'^, and require that the Kg daughter pions satisfy 
the same tracking criteria as the slow pions of the D*'^ 
analysis. The Kg decay vertex is required to be inside 
the beam pipe and to be well-separated from the D° de- 
cay vertex. These selection criteria yield an extremely 
clean Kg sample (> 99.5% pure), which is used to deter- 
mine fractional corrections to the overall magnetic field 
and to the energy losses in the beam pipe (E^^^) and, 
separately, in the SVT (E^^^^). The points represented 
as open squares in Fig. [T] show the strong dependence 
of the reconstructed Kg mass on laboratory momentum. 
Adjusting only the estimated energy losses and detector 
material flattens the distribution, but there is still a re- 
maining discrepancy. This discrepancy is shown by the 
open squares in Fig. [T] at high momentum and indicates 
an overall momentum scale problem. These two effects 
lead us to consider corrections to the laboratory momen- 
tum and energy of an individual track of the form 



TABLE I. Energy-loss and momentum correction parameters 
of Eq. @ which remove the momentum dependence of the 
reconstructed K% mass shown in Fig. [T] The nominal param- 
eters shift the average reconstructed masses to be the PDG 
mean value, also shown in Fig. [T] To estimate the associated 
systematic uncertainty, the procedure was repeated to give 
average reconstructed K^ masses itlcrpDG from the nominal 
value. 



Nominal 



For systematics 



'TIPDg(-K's) JTlPDG + IcTpDG mpDG — lo-pDG 



"bmp 
Osvt 



0.00030 
0.018 
0.059 



0.00030 
0.052 
0.059 



0.00019 
0.030 
0.059 



P 
E 






binp 



+ Obmp^ioss + Osvt^f 



(3) 



where the initial energy losses are determined by the 
Kalman filter based on the material model. We deter- 
mine the best set of correction parameters to minimize 
the x^ of the bin-by-bin mass difference between the 
7r"^7r~ invariant mass and the current value of the K% 
mass (mpDG {K%) ± Ictpdg = 497.614 ± 0.024 MeV) [ij. 

To assess the systematic uncertainty in values mea- 
sured from corrected distributions, we find new param- 
eter values by tuning the tt+tt" invariant mass to the 
nominal K^ mass shifted up and down by one standard 
deviation. These three sets of correction parameters are 
listed in Table HI The corrections corresponding to the 
nominal K^ mass remove the momentum dependence of 
m (tt+tt") as illustrated by the points shown as solid dots 
in Fig. m 

The best-fit value of a = 0.00030 corresponds to an in- 
crease of 4.5 Gauss on the central magnetic field. This is 
larger than the nominal 2 Gauss sensitivity of the mag- 
netic field mapping [ij]. However, the azimuthal depen- 
dence of Amo (discussed in Sec.p) indicates that the accu- 
racy of the mapping may be less than originally thought. 

The momentum dependence of Attiq in the initial re- 
sults is ascribed to underestimating the dE/dx loss in the 
beam pipe and SVT, which we correct using the factors 
&bmp (1-8%) and 6svt (5.9%). Typical dE/dx losses for 
a minimum ionizing particle with laboratory momentum 
2 GeV traversing the beam pipe and SVT at normal inci- 
dence are 4.4 MeV. The corrections are most significant 
for low-momentum tracks. However, the corrections are 
applied to all D*+ daughter tracks, not just to the slow 
pion. The momentum dependence was eliminated after 
the corrections were applied. All fits to data described in 
this analysis are performed using masses and lS.m values 
calculated using corrected 4-momenta. The MC tracks 



> 497.8 
" 497.7 



' 497.6 ^ 



497.5 
497.4 
497.3 
497.2 



497.1 



n Original • Corrected | | PDG K° mass ± 1 a 



■+++ 



++ 



^ 



f+ 



h^ 



{■4* 



jW 



,K'''\ 



*nn 



f 



10-^ 



1 

Lab momentum of slower tc [GeV] 



FIG. 1. Mass value of the K% obtained by fitting the invari- 
ant 7r'''7r~ mass distribution shown as a function of the slower 
pion laboratory momentum before and after energy-loss and 
momentum corrections have been applied. Note that the hor- 
izontal scale is logarithmic. 



are not corrected because the same field and material 
models used to propagate tracks are used during their 
reconstruction. 



FIT METHOD 

To measure T we fit the Ar/i peak (the signal) with 
a relativistic Breit-Wigner (RBW) function convolved 
with a resolution function based on a Geant4 MC sim- 
ulation of the detector response [18[. As in previous 
analyses [9[, we approximate the total D*^ decay width 
rTotai("^) ~ r£).D7r("T-), ignoring the electromagnetic 
contribution from £)*+ — >• D'^j. This approximation has 
a negligible effect on the measured values as it appears 
only in the denominator of Eq. ([T]). For the purpose 
of fitting the Am distribution we obtain dT{Am)/dAm 
from Eqs. ((T|) and ([2]) by making the substitution m = 
m{D'^) + Am, where m{D'^) is the current average mass 



of the D° meson [17 1 



Our fitting procedure involves two steps. In the first 
step we model the resolution due to track reconstruction 
by fitting the Am, distribution for correctly reconstructed 
MC events using a sum of three Gaussians and a function 
to describe the non-Gaussian component. The second 
step uses the resolution shape from the first step and 
convolves the Gaussian components with a relativistic 
Breit-Wigner of the form in Eq. ([1]) to fit the Am. dis- 
tribution in data, and thus measure F and Attiq. We fit 
the Am distribution in data and MC from the kinematic 
threshold to Am = 0.1665 GcV using a binned maximum 
likelihood fit and an interval width of 50 keV. 



Modeling experimental resolution 




0.155 



0.16 0.165 
A m [GeV] 



V ■••"..■ 


• ..• 


.^- 


"■ • v^i'',." ■. ' ^'V*'.. -■•■,', -0'-'.-^ 


■' , •' 




■%?•*.■ ./ -V-csJfi^'-' 


•■>^. 




0.14 



0.145 0.15 0.155 0.16 0.165 
A m [GeV] 



III 



"-V^A'--? : ^ 



.»» ■ ■ * — V « **.*' 



•.-^^.-^ 



FIG. 2. (color online) Binned maximum likelihood fit to the 
Am resolution distribution of MC samples for both D" decay 
modes. The interval size is 50 keV, and the high mass tails are 
dominated by low statistics. Normalized residuals are defined 
as (A^obaerved - A^'predicted) /-y/ ^predicted • The shapes iu the 
distribution of the normalized residuals are from dominance 
by Poisson statistics. In the peak region the total PDF is 
visually indistinguishable from the Gaussian component of 
the resolution function. 

We generate samples of Z3*+ decays with a line width 



of 0.1 keV, so that all of the observed spread is due to 
reconstruction effects. The non- Gaussian tails of the dis- 
tribution are from events in which the tt^ decays to a /i in 
flight and where coordinates from both the tt and fj, seg- 
ments are used in track reconstruction. Accounting for 
these non-Gaussian events greatly improves the quality 
of the fit to data near the Am peak. 

We fit the Am distribution of the MC events with the 
function 

fNGSNG{Am;q,a) + (1 - f^'§) [/iG(Am; //i,ai) 
+/2G (Am; /X2, 12) + (1 - /i - /2) G (Am; /is, a^)] 

(4) 

where the G {Am; ^i,ai) are Gaussian functions and 
/jvG) /i) /2 are the fractions allotted to the non-Gaussian 
component and the first and second Gaussian compo- 
nents, respectively. The function describing the non- 
Gaussian component of the distribution is 



Sng (Am,; q, a) — Am w' e° 



(5) 



where u = (A77i/A77ithios) — 1 and A?7ithrcs = 'Ti^r is the 
kinematic threshold for the D*^ -^ D'^tt^ process. For 
Am < Amthros, Sng is defined to be zero. 

Figure [2] shows the individual resolution function fits 
for the two Z?° decay modes. Each plot shows the to- 
tal resolution probability density function (PDF) as the 
solid curve, the sum of the Gaussian contributions is rep- 
resented by the dashed curve, and the Sng function as 
a dotted curve describing the events in the tails. The 
resolution functions should peak at the generated value, 
Am^^ = m(D*(2010)+)-m(LiO) ^. However, the fit- 
ted values of the /i; differ slightly from Amg^'-^ , with the 
result that the superposition of the three Gaussian func- 
tions reaches a maximum at a slightly larger value of Am. 
The Sng function is excluded from this calculation as the 
peak position is not well defined and Sng describes less 
than 1% of the signal. We take this reconstruction bias 
as an offset when measuring Attiq from data and denote 
this offset by Stuq. The SniQ offset is 4.3 keV and 2.8 keV 
for the D'^ -^ A'^7r+ and D'^ — > K^tt^tt^tt^ modes, re- 
spectively. As discussed in Sec. „ although the values of 
6m,o are larger than the final estimates of the systematic 
uncertainty for Amg, they are required for an unbiased 
result from fits to the validation signal MC samples. The 
systematic uncertainty associated with 6m,o is implicitly 
included when we vary the resolution shape, as discussed 
in Sec. „ 



Fit Results 

The parameters of the resolution function found in the 
previous step are used to create a convolved RBW PDF. 
In the fit to data, Sng has a fixed shape and relative 
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FIG. 3. (color online) The results of the fits to data for each 
D" decay mode. The fitted parameter values are summarized 
in Table [III The solid curve is the sum of the signal (dashed 
curve) and background (dotted curve) PDFs. The total PDF 
and signal component are visually indistinguishable in the 
peak region. 



fraction, and is not convolved with the RBW. The rel- 
ative contribution of Sng is small and the results from 
the fits to the validation signal MC samples are unbiased 
without convolving this term. We fit the data using the 
function, 



V{Am; e, F, Amo, c) — 
S{Am; e,r, Attiq) 



fs 



/5(Am)d(Am) 



+ (1 - fs) 



B{Am; c) 



jB{Am)d{Am) 



(6) 



where fs is the fraction of signal events, S is the signal 
function 



S{Am) = RBW® 
(l-/A^G)[/f^^G(Am;Mr 



Am, 



MC „MC 



+/f^G(Am;Mr 

+ (1- ./r^-/f^) 

G(Am;/i3^^^ 



A<^ 



hIC 



(1 



Aml'^^af^{X 






0) 



I fMCc ( \^. rMC ^MC\ 

+JNG ^NG(^rn,q ,a j, 
and B is the background function 



B{Am) — Am \fu 



u e 



(7) 



(8) 



where u = (ATTi/ATTithres) — 1. The nominal RBW func- 
tion pole position is located at tti = Atoq + m{D'^) and 
has natural line width T. The Gaussian resolution func- 
tions convolved with the RBW have centers offset from 
zero by small amounts determined from MC, fii — Atuq . 
The widths determined from MC, <jf^'^ , are scaled by 
(1 -|-e) where e is a common, empirically determined con- 
stant which accounts for possible differences between res- 
olutions in data and simulation. The validation of the fit 
procedure is discussed in Sec. n 

Figure [3] shows the fits to data for both D'^ decay 
modes. The total PDF is shown as the solid curve, 
the convolved RBW-Gaussian signal as the dashed curve, 
and the threshold background as the dotted curve. The 
normalized residuals show the good agreement between 
the data and the model. Table HIl summarizes the results 
of the fits to data for the two modes. The table also 
shows the fitted S/B at the peak and in the Am tail on 
the high side of the peak. The long non-Gaussian tail of 
the RBW is required for the model to fit the data so well. 



SYSTEMATIC UNCERTAINTIES 

Wc have measured the pole mass and the width of 
the D*'^ meson with unprecedented precision, analyz- 
ing a high-purity sample of continuum produced Z?*"*" 
in e"*"e~ collisions at approximately 10.6 GeV, equivalent 
to approximately 477 fb^^, collected by the BaBAR de- 
tector. The results for the two independent £)° decay 
modes agree with each other well. The dominant sys- 
tematic uncertainty on the RBW pole position comes 
from the azimuthal variation. For the decay mode D'^ — > 
K-Tr+ we obtain F = (83.4 ± 1.7 ± 1.5) keV and Amo = 
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TABLE II. Summary of the results from the fits to data for 
the D" — >■ K^n'^ and D'^ — ^ K~n^Tr~n^ channels (statistical 
uncertainties only); S/B is the ratio of the convolved signal 
PDF to the background PDF at the given value of Am, and 
f is the number of degrees of freedom. 



Parameter 


D" -^ Ktv 


D° -^ K-K-K-K 


Number of signal events 


138 539 ±109 


174 286 ±150 


r ( keV) 


83.5 ±1.7 


83.2 ±1.5 


scale factor, (1 + e) 


1.06 ±0.01 


1.08 ±0.01 


Amo ( keV) 


145 425.6 ±0.6 


145 426.6 ±0.5 


S/B at peak 
(Am = 0.14542 (GeV)) 


2700 


1130 


S/B at tail 
(Am = 0.1554 (GeV)) 


0.8 


0.3 


2 / 


574/535 


556/535 



tcrmincd from the corresponding MC subset. 

If the fit results from the disjoint subsets are compat- 
ible with a constant value, in the sense that x^ /v < 1, 
wc assign no systematic uncertainty. However, if we find 
X^ /v > 1 and do not determine an underlying model 
which might be used to correct the data, wc ascribe an 
uncertainty using a variation on the scale factor method 
used by the Particle Data Group (see the discussion of 
unconstrained averaging [17|]). The only sample which 
we do not fit to a constant is that for Attiq in intervals of 
azimuthal angle. We discuss below how we estimate the 
associated systematic uncertainty. 

In our version of this procedure, we determine a factor 
that scales the statistical uncertainty to the total un- 
certainty. The remaining uncertainty is ascribed to un- 
known detector issues and is used as a measure of sys- 
tematic uncertainty according to 



(145 425.6 ± 0.6 ± 1.8) kcV while for the decay mode 
D° -^ A'-7r+7r-7r+ wc obtain T = (83.2 ± 1.5 ± 2.6) kcV 
and Amo = (145 426.6 ± 0.5 ± 2.0) keV. Account- 
ing for correlations, we obtain the combined measure- 
ment values r = (83.3 ± 1.3 ± 1.4) kcV and Amo = 
(145 425.8 ±0.5 ±1.8) keV. 

We divide the data into disjoint subsets corresponding 
to intervals of laboratory momentum p of the D*+, az- 
imuthal angle cj) of the D*~^ in the laboratory frame, and 
reconstructed D" mass, to search for variations greater 
than those expected from statistical fluctuations. We also 
vary the resolution shapes according to their uncertain- 
ties and correlations, and estimate uncertainties related 
to our choice of end point for the fit, the nominal interval 
size, and the radius assumed in the Blatt-Weisskopf form 
factor. Finally, we assess systematic uncertainties asso- 
ciated with possible biases observed in fits to MC valida- 
tion samples. All these uncertainties are estimated inde- 
pendently for the £)" -^ A'^7r+ and D" -> Ar~7r+7r^7r+ 
modes, and are summarized in Table Hill 

After the corrections to the material model and mag- 
netic field, the laboratory momentum dependence of the 
RBW pole position is all but eliminated. We find that F 
does not display an azimuthal dependence, however Attiq 
does. Neither F nor Amo displays a systematic variation 
with reconstructed D° mass. We discuss in the next sub- 
section how we estimate systematic uncertainties using 
disjoint subsets. 



Systematics using disjoint subsets 

We chose to carefully study laboratory momentum, re- 
constructed D'^ mass, and azimuthal angle (p in order to 
search for variations larger than those expected from sta- 
tistical fluctuations. For each disjoint subset, we use the 
resolution function parameter values and Amo offset de- 



^y^ 



1 



(9) 



where the scale factor is defined as S^ = x^/^- The x^ 
statistic gives a measure of fluctuations, including those 
expected from statistics, and those from systematic ef- 
fects. Once we remove the uncertainty expected from 
statistical fluctuations, we associate what remains with 
a possible systematic uncertainty. 

We expect that x^/^ will have an average value of unity 
if there are no systematic uncertainties that distinguish 
one subset from another. If systematic deviations from 
one subset to another exist, then we expect that x^/^ 
will be greater than unity. Even if there are no system- 
atic variations from one disjoint subset to another, x^/^ 
will randomly fluctuate above 1 about half of the time. 
To be conservative, we assume that any observation of 
X^/t^ > 1 originates from a systematic variation from one 
disjoint subset to another. This approach has two weak- 
nesses. If used with a large number of subsets it could 
hide real systematic uncertainties. For example, if in- 
stead of 10 subsets we chose 1000 subsets, the larger sta- 
tistical uncertainties wash out any real systematic vari- 
ation. Also, if used with a large number of variables, 
about half the disjoint sets will have upward statistical 
fluctuations, even in the absence of any systematic vari- 
ation. We have chosen to use only three disjoint sets of 
events, and have divided each into 10 subsets to mitigate 
the effects of such problems. 

We chose the range for each subset to have approxi- 
mately equal statistical sensitivity. In Figs. 4(a) and |4(b)| 
we show the fit results in subsets of laboratory momen- 
tum for F and Amo, respectively. Neither D'^ mode dis- 
plays a systematic variation; however, we assign small 
uncertainties for each channel using Eq. ([9]). Similarly, 
Figs. 4(c) and |4(d)| show the results in ranges of re- 
constructed D° mass for F and Amo. While neither 
mode displays a systematic dependence, the width for 
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TABLE III. Summary of systematic uncertainties with correlation, p, between the D" — > K^iv'^ and Z?" — >■ K^tv^tv^tv^ 
modes. The K^tt^ and K^tv^tt^tv^ invariant masses are denoted by m (D^cco)- The methods used to calculate or define the 
correlations are described in Sec. □ The total systematic uncertainties are calculated according to the procedure defined in Sec. □ 



Source 


CTsya 


(r) [keV] 


P 


(Jsys 


(Amo)[keVl 


p 




K-K 


K-R-K-R 




K-K 


K-R-K-R 




p dependence 


0.88 


0.98 


0.47 


0.24 


0.20 


0.28 


m (-Drcco) dependence 


0.00 


1.53 


0.56 


0.04 


0.00 


0.22 


Azimuthal dependence 


0.62 


0.92 


-0.04 


1.65 


1.81 


0.84 


Magnetic field and material model 


0.29 


0.18 


0.98 


0.75 


0.81 


0.99 


Blatt-Weisskopf radius 


0.04 


0.04 


0.99 


0.00 


0.00 


1.00 


Variation of resolution shape parameters 


0.41 


0.37 


0.00 


0.17 


0.16 


0.00 


Am fit range 


0.83 


0.38 


-0.42 


0.08 


0.04 


0.35 


Background shape near threshold 


0.10 


0.33 


1.00 


0.00 


0.00 


0.00 


Interval width for fit 


0.00 


0.05 


0.99 


0.00 


0.00 


0.00 


Bias from validation 


0.00 


1.50 


0.00 


0.00 


0.00 


0.00 


Radiative effects 


0.25 


0.11 


0.00 


0.00 


0.00 


0.00 


Total 


1.5 


2.6 




1.8 


2.0 





mode is assigned its largest uncertainty 
of 1.53 keV using Eq. ®. 

In this analysis we observe sinusoidal variations in the 
mass values for D" -^ A'^7r+, Z?" — > A'^7r+7r^7r+, and 
K^ — >■ 7r+7r~. An extensive investigation revealed that 
this dependence most probably results from small devi- 
ations in the magnetic field from the map used in track 
reconstruction. For example, if during the field mapping 
(see Ref. ij] ) the propeller arm on which the probes were 
mounted flexed then the radial and angular components 
of the magnetic field would have mixed. The important 
aspect for this analysis is that, for such deviations, the 
average value is unbiased by the azimuthal variation as 
evidenced by the corresponding behavior of the K^ mass 
value. Figures 4(e) and |4(f)| show F and Amo, respec- 
tively, in subsets of azimuthal angle. The width does 
not display a dependence, but each mode is assigned a 
small uncertainty using Eq. ([9]). The lack of systematic 
variation of F with respect to is notable because Atoq 
(which uses reconstructed D masses) shows a clear de- 
pendence such that the results for the D^ — >■ K~'k'^ 
and _D" ^- A'^7r+7r^7r+ datasets are highly correlated, 
and shift together. The signs and phases of these vari- 
ations agree with those observed for D^ — > A'~7r+, 
D" ^ A'^7r+7r-7r+, and A'g -^ tt+tt". We take half 
of the amplitude obtained from a sinusoidal fit as a con- 
servative estimate of the uncertainty. 



Additional systematics 

We assess the uncertainty associated with the correc- 
tion parameters for the detector material model and mag- 
netic field by examining the variation between the nomi- 
nal parameter values and those obtained by tuning to the 
TTipDG (-^s) -'- I'^'PDG niass values |17| . The width mea- 



sured from the D^ — > A'^7r+ mode fiuctuates equally 
around the value from the fit using the nominal correc- 
tion parameters. We take the larger of the differences 
and assign an uncertainty of 0.29 keV. The value of Ajtiq 
for this mode fluctuates symmetrically around the nom- 
inal value and we assign an uncertainty of 0.75 keV. The 



width measured from the D 



R- 



fluctuates 



asymmetrically around the nominal value, and we use 
the larger difference to assign an uncertainty of 0.18 keV. 
The value of Amo for this mode fluctuates symmetrically 
around the nominal value, and we assign an uncertainty 
of O.SlkeV. 

We use the Blatt-Weisskopf radius r = 1.6 GeV^^ (~' 
0.3 fm) [13[. To estimate the systematic effect due to the 
choice of r we refit the distributions treating the Z?*+ as 
a point-like particle (r = 0). We see a small shift of F, 
that we take as the estimate of the uncertainty, and an 
effect on the RBW pole position that is a factor of 100 
smaller than the fit uncertainty, that we neglect. 

We determine the systematic uncertainty associated 
with the resolution function by refitting the data with 
variations of its parametrization. We take the covari- 
ance matrix from the fit to MC resolution samples for 
each mode and use it to generate 100 variations of these 
correlated Gaussian-distributed shape parameters. We 
use these generated values to refit the data, and take 
the root-mean-squared (RMS) deviation of the resulting 
fit values as a measure of systematic uncertainty. This 
process implicitly accounts for the uncertainty associated 
with the reconstruction offset. 

Our choice of fit range is somewhat arbitrary, and so 
we study the effect of systematically varying its endpoint. 
This changes the events associated to the RBW tails and 
those associated to the continuum background. This ef- 
fectively estimates the uncertainty in the background pa- 
rameterization. We expect that the larger the Am fit 
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FIG. 4. The values of F (left) and Amo (right) obtained ftom fits to data divided into 10 disjoint subsets in laboratory 
momentum p (top row), reconstructed D" mass (center row), and azimuthal angle (bottom row). The quantities p and 4> are 
defined by the D*'^ momentum. Each point represents an individual fit and each horizontal line is the nominal fit result (i.e. 
integrating over the variable). The correlation value of F (or Amo) measured from the D° — >■ K~tt'^ and D" — >■ K^tv^tt^tv^ 
samples for each of the variables chosen is given above each plot. The widths from the nominal fits and the weighted average 
agree well and the corresponding lines are visually indistinguishable. 
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range, the better description we will obtain for the back- 
ground and, thereby, a better description of the RBW 
tails. We vary the end point anyway and take 1 MeV 
steps from the nominal fit end point of Am = 0.1665 GeV 
down to Am = 0.16 GeV. Values below Am = 0.16 GeV 
arc too close to the signal region to provide a reasonable 
choice of end point. There is no clear way to estimate 
the associated systematic uncertainty, and so we take the 
largest deviation from the nominal fit as a conservative 
estimate of the uncertainty. 

The shape of the background function in Eq. ^ is 
nominally determined by only the parameter c and the 
residuals in Figs. 3(a) and |3(b)| show signs of curva- 
ture indicating possible systematic problems with the 
fits. Changing the end points over the range consid- 
ered changes the values of c substantially from —1.97 
to —3.57, and some fits remove all hints of curvature in 
the residuals plot. Wc also examine the influence of the 
background parametrization near threshold by varying 
the square root in Eq. ([8]) from 0.5 to nearby values 0.45 
and 0.55. The value of the fractional power controls the 
shape of the background between the signal peak and 
threshold. For example, at Am = 0.142 GeV changing 
the power from 0.5 to 0.45 and 0.55 varies the background 
function by -fl8% and -15%, respectively. The RBW 
pole position is unaffected by changing the background 
description near threshold while F shifts symmetrically 
around its nominal values. Wc estimate the uncertainty 
due to the description of the background function near 
threshold by taking the average difference to the nominal 
result. 

In the binned maximum likelihood fits we nominally 
choose an interval width of 50 keV. As a systematic check, 
the interval width was halved and the fits to the data 
were repeated. The measured F and Attiq values for both 
modes arc identical except for the width measured in 



the D 







K TT+TT 7r+ decay mode. Wc take the full 



difference as the systematic uncertainty for the choice of 
interval width. 



Validation samples 



identical between the fits to the total validation signal 
MC sample and the subsets, so we do not assign a sys- 
tematic uncertainty to the poles for radiative effects. The 
widths measured in each mode show a small difference to 
the results from the nominal validation sample. We take 
half of this difference as a conservative estimate of the 
systematic uncertainty associated with radiative effects. 



TABLE IV. Summary of results of the fits to the D° -^ K--K+ 
and D" — >■ K~Tv^'K~Tv'^ validation MC samples. The width 
from the D'^ — >■ K~n'^n~Tv^ decay mode has a small bias, 
which we take as a systematic uncertainty. 



Fit value 



Generated D° -^ Ktv D° -)■ Ktttvtt 



r[keV] 88.0 

scale factor, 1 + e 1.0 
Amo[keV] 145400.0 



88.5 ±0.8 89.5 ±0.6 

1.003 ± 0.004 1.000 ± 0.001 

145399.7 ± 0.4 145399.2 ± 0.4 

613/540 770/540 



Determining correlations 

The fourth and seventh columns in Table IIIII list 
the correlations between the D° — ?► K~tt'^ and Z?° — > 
K^tt^tt^t:^ systematic uncertainties. These correlations 
are required to use information from both measurements 
to compute the average. The correlations in laboratory 
momentum, reconstructed D^ mass, and azimuthal an- 
gle disjoint subsets are calculated by finding the cor- 
relation between the 10 subsets of Z?° — >■ A'~7r+ and 
_D" -^ Ar^7r+7r^7T+ for each of the variables. In a simi- 
lar way we can construct datascts using the sets of cor- 
rection parameters for magnetic field, detector material 
model, and the Am fit range. We assume no correla- 
tion for the resolution shape parameters and the vali- 
dation shifts, which are based on the individual recon- 
structions. Our studies show that the values chosen for 
the Blatt-Wcisskopf radius and interval width affect each 
mode identically, so we assume that they arc completely 
correlated. 



We generate signal MG with F = 88keV and Amp = 
0.1454 GeV. The background is taken from a MC cocktail 
and paired with the signal in the same ratio as from the 
corresponding fits to data. Fits to both decay modes 
describe the validation samples well. The fit results are 
summarized in Tabic IIVI We observe a small bias in the 
fitted width for the I?° -^ Ar~7r+7r~7r+ mode. To be 
conservative, we take the full difference between the fitted 
and generated value of the width and assign a 1.5 kcV 
error. 

We also investigated the uncertainty due to radiative 
effects by examining the subset of these events generated 
without PHOTOS fli|. The values of the RBW pole are 



Consistency checks 

In addition to the investigations into the sources of sys- 
tematic uncertainty, we also perform a number of consis- 
tency checks. These checks are not used to assess system- 
atics, nor are they included in the final measurements, 
but serve to reassure us that the experimental approach 
and fitting technique behave in reasonable ways. First, 
we lower the p* cut from 3.6 GeV to 2.4 GeV. This allows 
in more background and tracks with poorer resolution, 
but the statistics increase by a factor of three. Corre- 
spondingly, the signal-to-background ratios measured at 
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the peak and in the tails decrease by approximately a 
factor of three. The fit results for this larger dataset are 
consistent with the nominal fit results. The second con- 
sistency check widens the reconstructed D'^ mass window 
from 10 MeV to 30 MeV. Again, this increases the number 
of background events and improves statistical precision 
with central values that overlap with the nominal fit re- 
sults. Finally, we fix the scale factor in the fit to data 
to report statistical uncertainties on F similar to those 
in the measurement by CLEO [9| . Our reported "statis- 
tical" uncertainties on F are from a fit in which e floats. 
As expected, there is a strong negative correlation be- 
tween e and F with p (F, e) sa —0.85. If less of the spread 
in the data is allotted to the resolution function then it 
must be allotted to the RBW width, F. We refit the 
D^ — >■ K~Tr~^ and D'^ — )- K'tt^tt^tt'^ samples fixing e to 
the value from the fit where it was allowed to fioat. This 
effectively maintains the same global minimum while de- 
coupling the uncertainty in F from e. The statistical un- 
certainty on the width decreases from 1.7 keV to 0.9 keV 
for the D° — > K^tt^ decay mode and from 1.5 keV to 
0.8 keV for the I?° -^ K^tt^tt^tt^ decay mode. 



COMBINING RESULTS 

Using the correlations shown in Table IIIII and the for- 
malism briefiy outlined below, we determine the values 
for the combined measurement. For each quantity, F and 
Attiq, we have a measurement from the D'^ — > K^tt^ and 
D'^ — >■ K~Tr~^7T~TT~^ modes. So, we start with a 2 x 2 co- 
variance matrix 



V 



'K-n 



cov^Kir, Kmnr) 



2 _L 2 V^ 

•^KTT.stat + ^A'7r,sys l^i Pi ^Kir^i (^Kir-KTr^i 

E2 \ 2 

i pi O'Kn.i O'KTTTTTT.i ^K7T'!T7T,stat + '^i<:7r7r7r,sys 

(10) 

where i is an index which runs over the sources of system- 
atic uncertainty. In the final step we expand the nota- 
tion to explicitly show that the diagonal entries incorpo- 
rate the full systematic uncertainty and that the statis- 
tical uncertainty for the individual measurements plays 
a part in determining the weights. The covariance ma- 
trices are calculated using Table IIIII and the individual 
measurements. From the covariance matrix we extract 
the weights, w, for the best estimator of the mean and 
variance using Wi 



l^k ^ik I l^jk ^jk 



Wr 



WAmo — 



WKtt 



WKtt 
WKiTTT-n 



0.648 
0.351 

0.647 
0.353 



(11) 



(12) 



The weights show that the combined measurement is 
dominated by the cleaner D° — > K^tt'^ mode. The total 
uncertainty can be expressed as 

O-^ = ^ (WjCTstat.i)^ 
i=l,2 

i=l,2 i=l,ll 

(13) 

The statistical contribution is the first term and is 
simply calculated using the individual measurements 
and the weights. The remaining two terms repre- 
sent the systematic uncertainty, which is simply the 
remainder of the total uncertainty after the statisti- 
cal contribution has been subtracted. The weighted 
results are F = (83.3 ± 1.3 ± 1.4) kcV and Amo = 
(145 425.8 ±0.5 ±1.8) keV. 

SUMMARY AND CONCLUSIONS 

We have measured the pole mass and the width of 
the D*^ meson with unprecedented precision, analyz- 
ing a high-purity sample of continuum-produced D*^ 
in e^e~ collisions at approximately 10.6 GeV, equivalent 
to approximately 477 fb~^, collected by the BaBAR de- 
tector. The results for the two independent D^ decay 
modes agree with each other well. The dominant sys- 
tematic uncertainty on the RBW pole position comes 
from the azimuthal variation. For the decay mode D^ — > 
K-TT+ we obtain F = (83.4 ± 1.7 ± 1.5) keV and Amo == 
(145 425.6 ± 0.6 ± 1.8) keV while for the decay mode 
Z?" -^ K-1T+TT-TT+ we obtain F = (83.2 ± 1.5 ± 2.6) keV 
and Amo = (145 426.6 ± 0.5 ± 2.0) keV. Account- 
ing for correlations, we obtain the combined measure- 
ment values F = (83.3 ± 1.3 ± 1.4) kcV and Amo = 
(145 425.8 ±0.5 ±1.8) keV. 

The experimental value of go'D-n is calculated using 
the relationship between the width and the coupling con- 
stant. 



F (i?%+) 
F (i5V+) 



F (i?+^") 
F (Z?+^0) 



F i?+7 



9d*D0t:+ 3 I ffcD+Tr" 3 

247rm2 ./'^^ ^ 24^m?, ^^" 



(14) 
(15) 

(16) 



'£>* + 



'_D' + 



where we have again ignored the electromagnetic con- 
tribution. The strong couplings can be related through 
isospin by gD*D°-K+ = —V^qn-'n+irO (9[. Using F and the 
mass values from Ref. [17| we determine the experimen- 
tal coupling g^?DO^+ = 16.92 ± 0.13 ± 0.14. The univer- 
sal coupling is directly related to the strong coupling by 
g = gD'D°7T+ Itt/ (2-^m£)m£)») . This paramctrization is 
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TABLE V. Selected rows from Table 11 of Ref. [21(. State 
names correspond to the current PDG listings. The third 
column is the ratio, R = T/cf, extracted from the model 
in Ref. [2l||. The values of g were obtained from the data 
available in 2001. 



State 



Width (r) 



R 

(model) 



D* (2010) + 
Di (2420)° 
D2 (2460)" 



96±4±22keV 143 keV 0.82 ±0.09 



18.9 



+4.6 



MeV 



23 ± 5 MeV 



16 MeV 1.09: 

38 MeV 0.77 ±0.08 



1+0.12 
-0.11 



TABLE VL Updated coupling constant values using the latest 
width measurements. Ratio values are taken from Table [V] 
Significant differences are seen among the coupling constants 
calculated using the updated width measurements. 



State 


Width (r) 


R 

(model) 


9 


D*(2010) + 
Di (2420)" 
Dl (2460)" 


83.3 ± 1.3 ± 1.4 keV 

31.4 ±0.5 ±1.3 MeV 

50.5 ±0.6 ±0.7 MeV 


143 keV 
16 MeV 
38 MeV 


0.76 ±0.01 
1.40 ± 0.03 
1.15 ±0.01 



different from that of Ref. [9[ and is chosen to match a 
common choice when using chiral perturbation theory, as 
in Refs. ^^. With this relation and ,U = 130.41 MeV, 
we find g'^^'P = 0.570 ± 0.004 ± 0.005. 

The paper by Di Pierro and Eichten 2l| quotes re- 
suhs in terms of a ratio, R = T/g'^, which involves the 
width of the particular state and provides a straightfor- 
ward method for calculating the corresponding value of 
the universal coupling constant within their model. The 
coupling constant should then take the same value for 
the selected D^*^ decay channels listed in TableEl which 
shows the values of the ratio R extracted from the model 
and the experimental values for F, as they were in 2001. 
At the time of publication, g was consistent for all of the 



modes in Ref. [21|. In 2010, BaBAR published much more 
precise results for the Di (2420)° and D^ (2460)° ^. 
Using those results, this measurement of F, and the ratios 
from Table |Vl we calculate new values for the coupling 
constant g. Table IVII shows the updated results. We 
estimate the uncertainty on the coupling constant value 
assuming or ^ F. The updated widths reveal significant 
differences among the extracted values of g. 

The order of magnitude increase in precision confirms 
the observed inconsistency between the measured £)*+ 
width and the chiral quark model calculation by Di 
Pierro and Eichten 2l|. The precise measurements of 



the widths presented in Table IVII provide solid anchor 
points for future calculations. 
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